REPORT No. 888 


A THEORY OF UNSTAGGERED AIRFOIL CASCADES IN COMPRESSIBLE FLOW 

By Robert A. Spuer and H. Julian Allen 


SUMMARY 

By use of the methods of thin airfoil theory, which include 
effects of compressibility, relations are developed which permit 
the rapid determination of the pressure distribution over an 
unstaggered cascade of airfoils of a given profile, and the deter- 
mination of the profile shape necessary to yield a given pressure 
distribution for small churdlgap ratios. For incompressible 
flow the results of the theory are compared with available 
examples obtained by the more exact method of conformal trans- 
formation. Although the theory is developed for small chordj 
gap ratios, these comparisons show that it may be extended to 
chordigap ratios of order unity, at least for low-speed flows. 
Choking of cascades, a phenomenon of particular importance in 
compressor design, is considered. 

INTRODUCTION 

The wider use of gas turbines and other devices employing 
axial-flow compressors has increased the need for com- 
pressors with a liigh pressure rise per stage. In order to 
achieve this purpose it is necessary to use high velocity 
flows, thus increasing the possibility of losses through com- 
pression shock. A method is therefore desirable which wUl 
permit the design of compressor blades which have high 
critical compressibilitj" speeds. This result can be accom- 
plished if a cascade of airfoils representing the flow can be de- 
signed to give a desirable airfoil-section pressure distribution. 

This report attacks a portion of the problem by finding the 
relation between the profile shape and pressure distribution 
over airfoils arranged in an imstaggered cascade through the 
use of the approximate methods of thin airfoil theory origi- 
nally presented by Glauert in reference 1 and further de- 
veloped by the NACA in reference 2. The flow over an 
airfoil in cascade is related to that over a single airfoil in a 
free stream. The problem of finding the pressure distribu- 
tion over an airfoil in cascade or the shape of an ahfoil in 
cascade to give a required pressure distribution then reduces 
to the analogous problem for a single airfoil, which can be 
solved by known methods. 

The analysis involves the assumption that the gap between 
airfoils is large compared to the chord length. In particiflar, 
expressions relating the characteristics of a cascade airfoil to 
those of a free airfoil are expanded in a power series in cfg, 
where c is the chord and g the perpendicular chstance between 
ahfoil chord lines in the cascade, and powers of e/p higher 


than the second are neglected. Definitions of the sjonbols 
used are foimd in Appendix A. 

THEORY 

Consider an infinite unstaggered cascade of identical two- 
dimensional airfoils, as represented in figure 1 . The config- 
uration is specified by the chord/gap ratio cjg, the camber- 
line shape and thickness chstribution of the individual 
airfoils. For an individual ahfoil, the incident velocity is 
Uwith a corresponding density of p and angle of attack a'. 



The pressure distribution on a typical cascade ahfoil will 
be compared with that which would be obtained over a 
single ahfoil of the same shape in a stream of velocity Y and 
density p. The angle of attack a of the single ahfoil will not, 
in general, be the same as the angle of attack a of the 
ahfoil in cascade, and an e.xpression relating the two angles 
win be given. In the analysis to follow, primes will be used 
to designate properties of the cascade ahfoil. 

AERODYNAMIC CH.AEACTERISTICS OF A CASCADE AIRFOIL OF A GIVEN 

PROFILE SHAPE 


Effect of camber. — The cbordwise lift distribution for the 
free ahfoil is given by the Kutta-Joukowski relation 


dL_ „dV 
dx ^ dx 


( 1 ) 


5ol 
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where L is lift per unit span, T is the vorticity per unit span, 
and X is the chordwise distance from the leading edge to the 
to the point in question. The lift-distribution coefficient P 
is then defined as 


p=i^=-? ^ 
dx V dx 


( 2 ) 


In Appendix B, it is .shown that the stream velocity at 
the cascade airfoil is uniform so that the lift distribution 
for the cascade airfoil may be expressed as 


dx 


(3) 


and the lift-distribution coefficient of the cascade airfoil 
referred to the dynamic pressure <7 is 


g Wx V dx 


From equations (2) and (4), 


^ ‘"^~V\dx lx 


(4) 


(5) 


It is now desirable to express equation (5) in a more con- 
venient form. It is assumed that the vorticity distributions 
of the free and cascade airfoils, respectively, maj' be rep- 
resented by the following series (references 1 and 2): 


dV / 1 ® 

^= 2 F cot 2 sin 


n6 


dV 


-^-■2V cot ^d+^A„' sin nd^\ 


( 6 ) 


The new parameter 6 is related to x by the equation 
a; =2 (1 — cos 6) 

Equation (5) then becomes 


(7) 


P*— P=4 ^(Ao' -Ao) cot I e+'^iAl—An) sin nflj (8) 

The coefficients in equation (8) will now be evaluated by 
considering the conditions of flow at the airfoil boundaries. 
Let V and 1 be the vertical components of velocity induced, 
respectively, by the free airfoil and a particular cascade 
airfoil, and let Av be the vertical component of velocity 
induced by the other airfoils of the cascade. For small 
angles of attack, in order that the flow be tangent at the 
surfaces, the slope dyjdx of the airfoil camber line (which is 
the same for both free and cascade airfoils) miist be given 
by the following relations: 


dye ,v 

, v'Av 

+V+TJ 


(9) 


Appenclix B gives the vertical components of velocity 
induced by the vorticity and source-sink distributions which 
represent the airfoils of a cascade. For simplicity, the 
foUowmg symbols are used; 


48 

X=Vl— 


(Bll) 

(B4) 


where AI is the Jvlach number of the undisturbed stream. 
The velocity components calculated in Appendix B are 


II 

/ <x> \ 

— Ao+SA„ cos nd ) 

(B5) 

V n=l / 


II 

—A/+^,An cos ) 

n=l / 

(B6) 


and neglecting terms involving <t to the second and higher 
powders 

[(Ao'+|a/)-(2Ao'+AiO cos e] (B25) 

Substituting these relations in equations (9), there is 
obtained 

X (^-Ao+SA„ cos 

•^^=a'-|-X Ao'-h2jA„' cos ^ (10) 

-2|[(Ao'+^A/)-(2Ao'+A/) cos e] 

Since equations (IQ) are equal trigonometric series, the 
coefficients of cos nd can be equated: 

«-XAo=:.'-XA/-2| (Ao'+|a/)1 


( 11 ) 


XAi = XA/-f-2 (2 Ao'-|-A/) 

A2 = A/ 

A„=Ak (n?^l) 

From the first of equations (11), it is seen that the quan- 
tity (Ao' — Ao) depends on the angle of attack a of the free 
airfoil. This angle may be arbitrarily defined so that 


Ao'-Ac=0 


( 12 ) 


Combining equations (11) and (12), the following relations 
are obtained: 

A,/-Ao=0 


Ai' — Ai 2 —2 (2 Ao'-1-Ai') 

Ao' — A2=0 


and 


A/-A„=0 {ri9^l) 
a'-a=2j(Ao'+jA2') 


(13) 


(14) 
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]Making use of equations (13), equation (14) may also be 
written as follows: 


a'-a=2^(.4o+| a) 
From equations ( 8 ) and (13), there is obtained 


P*-P=-S ^ ( 2 ^o' + -d/) sin e 

A." 


(15) 


(16) 


Changing to the unprimed coefBcients with the use of equa- 
tions (13), and neglecting terms in a^, the expression becomes 


P*—P=—8^ ( 2 J.Q-i-x-li) sin d 


(17) 


It should not be interpreted from equations (16) and (17) 
that Ai' and are equal. Equation (17) is the result of 
dropping terms in cr and not of equating A/ and The 
difference between A/ and to the order cr is still given 
by equations (13). 

The section lift coefficient for the free airfoil is 




(18) 


Substituting from equations ( 2 ), ( 6 ), and (7) and performing 
the integration, one obtains 

Ci=t{2A(}-\-Ai) (19) 

Similarly, the quarter-chord moment coefficient is 


From equations (19) and ( 20 ) it follows that 

Cl 


( 20 ) 


2.do+^i = 


0 i 2 V 


Ci-r4c, 




( 21 ) 


With use of equations ( 21 ), equations (15) and (17) then 
become 


d — a=^( Cj-i-4c 




and 


drX W' 

8 <rCi 


■k\^ 

l&crCi 


sin 6 




K'-f) 


( 22 ) 

(23) 

(24) 


Effect of thickness. — From Appendix B, the horizontal 
velocity at an airfoil in cascade is greater than that at a 
single airfoil by the amount Au, which is given by the relation 


Au_Acr 


(B24) 


The quantity A is a function of airfoil thickness given by 
equation (B 20 ) and tabulated for various airfoils in table I, 
which has been taken from reference 3. 

The increment of velocity given by equation (B24) can be 
added to the velocity of the midisturbed stream Tq to give the 
true incident velocity Th- 


V=Vi-\-Au 

From equations (B24) and (25), 
Vi_ 


(25) 


(26) 


The density pi far ahead of the cascade is given by the 
following equation, derived for isentropic flow: 




-[(W 




1 


where 7 is the ratio of specific heat at constant pressure to 
that at constant volume. Using the binomial e.xpansion and 
neglecting terms containing powers of [{VilV}- — 1 ] higher 
than the first, there is obtained to the order of a 


Pi 


= 1 - 


A<rAP 

X® 


(27) 


The ratio of the dynamic pressure g at the center of pres- 
sure to the dynamic pressure gi far ahead of the airfoil is then 
found to be 

|=l-h 2 ^ (28) 

where 


1 = 1 -^ 


(29) 


The hft-distribution coefficient P' referred to the dynamic 
pressure gi is given by 


p/—i p* 

2 i 

The cascade section lift coefficient is 




(30) 


( 31 ) 


Substituting equations (24) and (28) in equation (31) and 
performing the integration, there is obtained, neglecting 
terms in <r, 

(32) 


r' r--2^ 
Cl — C ; — — ^ 




Equation (32) gives the relation between the lift coefficient 
of a cascade airfoil at angle of attack a' and that of a single 
airfoO. at angle of attack a. The relation between o' and a is 
given by equation ( 22 ) . 

The Mach number AP must now be related to the Mach 
number AI for the corresponding single airfoil. The velocity 
of sound is related to a by 




- 

“V T' 


( 33 ) 
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For isenfcropic flow 




which becomes, using equation (26) 


F,=r[i+(7-i)M^|^ 


Hence equation (33) becomes, to the order of 


and so 


V K / tti 


Using equations (26) and (35) then 


where 




.=l+(3^). 


Determination of pressure coefficient. — The. following co- 
efficients are useful in expressing the pressure on. an airfoil 
surface : 

Pt='^^ ■ ( 38 ) 

where pi is the local static pressure on the surface of the air- 
foil, and H, p, and q are, respectively, the total head, static 
pressure, and d3mamic pressure of the undisturbed stream. 
The variation of H with Mach number, assuming that the 
ratio of specific heats is 1.4, is given (reference 4) by the 
equations 

//=p-l-g(l + ,j) \ 

14 . ? ( 39 ) 


From these equations 


Si=\ + r}—P t 


A graph of as a function of Mach number is given in 
figure 2. The subscript I in equations (38) and (40) vfill be 
replaced bj^ U, L, or / to denote, respectively, the upper 
aud lower surfaces of an airfoil and the surface of its sym- 
metrical base profile. 

According to reference 5, the velocitj^ V'f along the base 
profile of a single airfoil is given for mcompressible flow by 
the expression 

^/= (41) 


where the pressure coefficients are referred to q. That is, 
p Pu—P - p PL-P .p' Pf-P 


and p is the static pressure of the stream corresponding to q. 



O .2 .4 .6 .8 /.O 

Atoc/j number, M 

Figure 2.— CompressibHiiy factor, as a functioa of Mach mimber. 


The velocity V/ along the base profile of the airfoil in a 
cascade in terms of T’'is siraplj’- 

V/_V; 

so that 

(42) 

since the cascade base-profile pressure coefficient P* is 
referred to the dynamic pressure q and to the static pressure 
p rather than to the corresponding quantities in the un- 
disturbed stream. Analysis (reference 3) shows that equa- 
tion (42) holds also for compressible flow, to an approxima- 
tion of the same order of magnitude as others already made. 

The upper- and lower-surface pressure coefficients referred 
to q are from reference 5. 

Pv*=l- 


Pi*=l- 

Althougli equations (43) were derived for incompressible 
flow, it can be shown that they are applicable in the case of 
flow in a compressible stream. 
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The coefBcients Pu* and Pl* can be converted into the 
coefficients Pu' and Pl referred to qi (the dynamic pressure 
in the undisturbed stream) and the corresponding static 
pressure by means of the foUovring equations, which are 
based on equations (28) and (40) : 


Reynolds number should be used but a theoretical distribu- 
tion (e. g., by references 2 , 6 , 7, or 8 ) modified by the Glauert- 
Prandtl or Karman-Tsien rule (references 9 and 10 ) is 
satisfactory. 

4. P* is found from equation (24). 


ScP=l-l-7}-Pu* 

SL^=l-hv-Pz.* 

PP=l-i-v'-SL' 

when 5 j' is the value corresponding to M' given bv equation 
(36). 

DETEE.MIN.4TIOX OF PRESSURE DISTEIBUTIO.N' FOR A GIVEN' PROFILE 
SH.APE IN CASC.ADE 

A method can now be outlined for finding the pressure 
distribution over an airfoil of given profile in an unstaggered 
cascade. It will be assumed at the outset that the Mach 
number M' of the approaching air is known and that the 
desired hft coefficient c/ of each of the cascade airfoils is 
given. 

1 . The Mach number of the corresponding free airfoil, M, 
is determined from equation (36) to a first order in <r as 


P^=P 


16gC; f ^ 

V c I, c 


5. is calculated from the free airfoil pressure 

distribution from equations (41) and (42) 


Vl-P/=Vl-P,= 


\T — PwT Vl — Pi 


(49) 


6 . ( 1 — P/*) and P* are combined to give ( 1 — Py*) and 
(1— Pi*) from equations (43) 


1 -Py* 

1-Pi* 


(l-P/ + iP*}^ 
' 1 — P/ 


(1 — P/* — -^P*)’^ 


(43) 


7. The cascade pressure coefficients Py' and Pi' are then 
found by means of equations (44): 

^u*=l-Pi7* + ^ 

Si* = l-Pi*-f, 


(45) 

wherein 

■’-SGJ 

X= Vl-(-4/')' (46) 

(Mr (47) 


and A may be obtained from table I or from equation (B20). 

2 . The lift coefficient of the corresponding free airfoil, 
Cl, is found from equation (32) by neglecting second-order 
terms in g as 


Ct=C 


rTl 




■')] 


(48) 


wherein ^ is given by equation (29). 

For the given airfoil the quarter-chord pitching moment 
and the angle of attack a corresponding to the lift 
coefficient Cj will be known and the angle of attack for the 
airfoil in cascade a (degrees) can be found from 


a' = a -L 57.3 ^ (Ci+^C„J (22) 

3. P—Pl—Pu for compressible flow at the Alach number 
AI is found for the single airfoil at the hft coefficient Ci. 
Preferably an e.vperimental pressure distribution at the 
appropriate ilach nuniber and at approximately the same 



(44) 


where t? and 17 ' corresponding to M and A/', respectively, are 
found from figure 2 . 


DETEE.MINATIO.N OF PROFILE SHAPE FRO.M A GIVEN PRESSURE 
DISTRIBUTION IN CASCADE 

The procedm'e for obtaining the profile shape for a given 
pressure distribution in cascade involves a process which 
is essentially the reverse of that just outlined. First the 
required pressure distribution is drawn. The initial choice 
must be made skillfully with reliance on experimental pres- 
sure distributions so that the distribution chosen shall 
correspond or nearly correspond to that obtainable with a 
real airfoil. From the chosen distribution the coefficients 
c{ and are determined by graphical or numerical 

integration and the pressure coefficients Pu' and Pp are 
read off at selected points. The method then consists of 
the following steps: 

1 . The Alach number M is determined from 

( 45 ) 
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2. Pu* and Pl* are found by applying equations (44) in 
reverse order, neglecting terms involving a- to the second and 
higher powers. 

Su'=l + v'-Pu 
Sp = l —Pl,' 

/ . r 

P + — Sxj'' 

Pp=l-\-yi—Si,* 


where r\' and ri correspond to M' and M. (See fig. 2.) 

3. (1— P/) is found by means of the following equations, 
which are readily obtained from equations (49): 


Vl-P/=Vl-P/= 


Vl-Pt;* + Vl-PP^ 


(49) 


The base profile should now be checked to see if it satisfies 
the closing condition given in reference 2, If the assumed 
pressure distribution does not correspond to a closed shape, 
it must be modified so that it docs. (See reference 2.) 

4. The single-airfoil lift-distribution coefficient P is found 
from the relations 


P* = Pi,^-Pu* 

(51) 


(52) 


.5. The final upper- and lower-pressure coefficients are 
calculated bj" means of the following equations from equa- 
tion (43); 



6. The angle of attack of the single airfoil is then given 
in degrees by 

a=a'-57.3 — (c/-i-4c,„^,/) (54) 

which is a modification of equation (22). 

7. The problem of finding a profile shape which will have 
a given pressure distribution in cascade has now been re- 
duced to the analogous problem for a single airfoil. It is 
now possible to calculate the corresponding incompressible- 
flow pressure distribution by the method of reference 10, or 
simply by multiplying Pxj and P^ by -yll—M-. An airfoil 
shape can be designed to give this pressure distribution by 
the method of reference 2. In the application of this 
method, it is generally necessary to make some changes in 
the pressure distribution, as previously noted, in order to 
make it correspond to a possible distribution for an actual 
airfoil. If the initial distribution has been well chosen, 
however, the changes will be minor. 


THE CHOKING OF CASCADES 

In the compressible adiabatic flow of a fluid in an elemen- 
tary stream tube of varying area Ai, the mass flow must be 
constant so that the logarithmic derivative 

dpi.dVi dAi 
p, IT + A, 

must vanish. Now the density pi and the velocity T', are 
related to the pressure pi by Bernoulli’s equation. 

^=^X~=-V,dVi 

Pl apt Pl 


The quantity dpijdpi is, of course, the square of the velocity 
of sound so that Bernoulli’s equation may be written 

dpi VdVi ^,,dVi 
pi~ a? - ^ 


where Mi is the local Mach number. Combining this 
expression with that of the logarithmic derivative then 


(1-xM;-) 


dAi 
A I 


From this relation it is seen that for a subsonic flow the 
area must decrease for a velocity increase, while at supersonic 
speeds the area must increase for a velocity increase. I’llien 
the Mach number is unity, then dA=0, so that sound speed 
is only attained w-here the area is a minimum. 

Considering the flow through the cascade of figure 1 as 
essentially unidimensional, then it is apparent that if the 
flow past the plane hh attains sonic speed, the mass flow 
through the cascade cannot be further increased and the 
cascade flow may be said to be “choked.” Of course, the 
flow tlirough a cascade is not unidimensional, but experi- 
ence wdth the similar phenomena of choking in wind tumuds 
has mdicated that the assumption of unidimensionality of 
flow' yields calculated choking Mach numbers in good agree- 
ment with experiment. 

For a unidimensional flow' it is shown in reference 3 lliat 
the ratio of the free area of the undisturbed stream A to the 
minimum floiv area Am is related to the choking ivlach 
number of the free stream Mch by 


A 1 

Am Mch 


1 



"12(7-1) ■ 


■ 1 ) 


(55) 


For the approaching stream wdiicli w'ill pass between any two 
airfoils of the cascades the free cross-sectional area is 


A=g cos a' 


w'hile the minimum area between the two airfoils is 


and hence 


Am=g—t=g 



Am 


cos a' 



where # is the maximum thickness of an airfoil. 


( 56 ) 
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Properly the thickness of the boundary layers on each 
surface of the airfoil should be added to the geometric 
thickness to obtain an effective thickness Using this 
value with equations (55) and (56) and a value of y (for air) 
of 1.4 then 


tj 

9 


= 1 


Mch cos a' 



(57) 


A practical estimate of the choking klach number may be 
found by assuming the angle of attack to be so small that 
cos ^ 1 and that the effective thickness is the geometric 
thickness. With these assumptions, values of the thickness/ 
gap ratio as a function of the choking Mach number are 
given in figure 3. 



Figube 3.— Choking Mach number ss a function of the ratio of airfoil thickness to 

cascade gap. 


It is shown in reference 3 that there exists another pos- 
sibility that choking may occur in the w'ake of the airfoils 
for very thin profiles as the result of the action of viscosity. 
In most cases this type of choking wiU not be of practical 
importance. 


DISCUSSION 

In order to check the accuracy of the equations which have 
been developed, pressure distributions for incompressible 
flow have been calculated at three lift coefficients for an un- 
staggered cascade of NACA 4412 airfoils. These cases were 
chosen to permit comparison with those which have been 
determined in reference 11 by the method of conformal 
transformation. The chord/gap ratio in all cases is 1.03 
and the hft coefficients considered are 0, 0.5, and 1.0. This 
comparison subjects the approximate theory of this report 
to a rather severe test, since the analysis has been developed 
on the basis of a small chord/gap ratio. The comparison 
with reference 11 indicates the agreement is good as may be 
seen in figures 4, 5, and 6. The single-airfoil pressure 
distributions, as obtained by conformal transformation in 
reference 11, are shown in the figures for the same lift 
coefficients as for the cascade airfoils. (In Appendix C 
the calculations to obtain the pressure distribution of figure 



Chordwise station, percent c 


C 

Ftgitee 4. — Pressure distribution for NACA 4412 airfoil alone and in cascade for "=l-03 and 

cj'=0. 
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Figure 5.— Pressure distrubution for NACA 1412 aiifoil alone and in cascade for -=1.03 and 


6 are given in detail. Table II gives all the necessaiy com- 
putations and serves to demonstrate the simplicity of the 
method.) It is evident from the figures that one effect of 
cascading airfoils is to impose a “negative camber in- 
fluence” on the pressure distributions as the lift coefficient is 
increased. Accordingly, for airfoils for use in cascade the 
camber must be exaggerated if a certain desired camber 
effect on the pressure distribution is to- be obtained. It is 
of interest to note that the calculated angles of attack are 
also in reasonably good agreement as ma}’- be seen in the 
following tabulation: 


Cl' 

Airfoil in cascade 

Airfoil alone 

a' (from ref- 
erence 11) 

a' (method of 
report) 

cc 

0 

-5. 9° 

' -6. 0° 

-4. 3° 

. 5 

-f I. 8° 

-f2. 3° 

-. 1° 

1. 0 

-1-9. 7° 

-h 10. 4° 

4-4. 0° 



Figure 6.— Pressure distribution for NAC A 4412 airfoil alone and In cascade fur ^ L03 and 

ci'=1.0 

Unfortunately no pressure distributions over unstaggered 
cascades at liigh airspeeds could be found so that the valiflity 
of the compressibility corrections developed in this report 
could not be determined. They have been developed, how- 
ever, in much the same way as the compressibility corrections 
of reference 3, which have been found to bo in good agree- 
ment with experiment. It is clear, however, that as the 
Mach number is increased the accuracy of the calculations 
wiU diminish unless the chord/gap ratio is simultaneously 
decreased. This will be particularly true as the choking 
Mach number is approached. 

One matter of interest concerns the effect of compressi- 
bihty on the lift-curve slope of a cascade of NACA 4412 
airfoils having chord/gap ratios of 0 (i. e., lone airfoil), 0.5, 
and 1,0 which is shown in figure 7. It is seen that the 
importance of the interference effects of the airfoils in cas- 
cade increase so rapidly with Mach number that, contrary 
to the usual expectation, a decrease in lift-curve slope with 
Mach number is indicated for high solidities. 

An examination of the equations for the choking of cas- 
cades indicates that care must be e.xercised with an unstag- 
gered cascade to keep the airfoil thickness small for high 




Lifi-curve slope ratio. 
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solidities if relatively Mgli free-stream Mach numbers are 
employed. Even for the cascade considered in figures 4, 5, 
and 6, for example, the choking Mach number as obtained 
from figure 3 is only about 0.64 based on the geometric 
thickness of the ahfoils. 

The existence of a boundary layer on the airfoil surfaces 
would, of course, increase the effective thickness which would 
reduce the choking Mach number. TVlien the pressure 
gradients are strongly adverse, as in the case of the airfoil 
of figure 6, the boundary-layer growth will be increased and 
its effect on the choking Mach number will be more pro- 
noimced. The proper choice of camber consistent with the 
design lift wiU serve to reduce the sharp pressure peak, 
thereb 3 ^ improving the critical as well as the choking Mach 
number. 

Aims Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., September, 1947. 



APPENDIX A 


LIST OF SYMBOLS 

The following symbols are used through this report: 
a velocity of sound 

A„ Fourier coefficients (See equation (7).) 

Bn Fourier coefficients (See equation (Bl4).) 

c airfoil chord 

Cl section lift coefficient 

section quarter-chord-moment coefficient 
g distance perpendicular to chord between airfoils in 
cascade 
H total head 

L lift per unit span 

M Mach number 

Mch choking Mach number 
p static pressure 

P without subscript, local lift at any chord station in 
coefficient form; vutb subscript, local pressure co- 
efficient (See equation (38).) 
q dynamic pressure 

Q source, strength 

r radial distance in polar coordinates 
S local pressure, coefficient (See equation (38).) 
t maximum airfoil thickness 

T absolute temperature 

u horizontal component of velocity 

V vertical component of velocity 

A» total induced vertical velocity at the airfoil under 
consideration due to the other airfoils in the cas- 
cade 

V velocity 

X coordinate of points on chord line as measured from 
leading edge 

Pc ordinate of mean camberjine 

Pt ordinate of base profile 

a. angle of attack 

S60 


7 ratio of specific heat at constant pressure to specific 
heat at constant volume (Cp/ed 
r circulation per unit span 

/ 

$ compressibLiity correction factor ( 1 — g- ) 

i\ compressibility factor (See equation (39) and fig, 2.) 

6 angular coordinate of points on chord line (Sec 

equation (7).) 

X compressibility correction factor (Vl—M*) 

A factor depending upon shape of base profile (See 

equation (B20) and table I.) 

/I compressibility correction factor ^1 ^ ^ 

p mass density 

<r factor depending upon solidity of cascade ( 

r inclination of the radius r to a line normal to the air- 
foil chord (See fig. 1.) 

# . polar angle in polar coordinates (positive counter- 

clockwise) 

SLPERSCKirrS 

' used to distinguish properties of cascade airfoil 

* denotes cascade airfoil characteristics as coefficients 

referred to dynamic pressure of incident stream at 
center of pressure 

SUBSCRIPTS 

‘ denotes values in stream far ahead of cascade (e.vccpt 
when used to indicate a numbered Fourier coefficient) 
I denotes local conditions at point in fluid 
L denotes values on lower surface of airfoil 

U denotes values on upper surface of airfoil 

j refers to base profile 
V denotes a velocity induced by a vortex 

s denotes a velocity induced by a source or sink 



APPENDIX B 


VELOCITY COMPONENTS INDUCED BY AIRFOILS IN CASCADE 

For incompressible flow over a single airfoil, the vertical 
velocity (perpendicular to the chord) induced at the point on 
the chord Xa by its own vorticity distribution is 


2vJo x—Xo 


(Bl) 


To obtain the velocity for compressible flow, the factor 


1— sin^ $ 


(B2) 


must be apphed to the right-hand side of equation (Bl) as 
is shown in reference 9. Here A/ is the Mach number of the 
undisturbed stream and $ is the angle between the stream 
direction and a line drawn from the vortex to the point in 
question. For a single airfoil is always close to 0° or 180° 
and the induced velocity in a compressible stream is ap- 
pro.ximately 


A 


/l-AP 


where 


dr , 

pdr“ X 1 

dr , 
p'dr 

(B3) 

Jo x—Xa 2xJ 

^0 X — Xo 



(B4) 


Equation (B3) can be integrated by substituting for dVjdx 
and for x the following: 

^=2F^*4q cot sin n6^ (6) 

(1— cosS) (7) 


a line from the vortex to the point in question, as indicated 
in figure 1. The factor in parentheses is that necessary to 
give the correct result for compressible flow. It is derived 
from the factor (B2) with use of the relation sin #= — cos r. 
Since r^={mg(cos rj and tan Tr^=ix—x^)lmg so that 


sm T jTj 


and 


cos 


X — Xg 

mg 



the preceding equation may be written 




dV jdx 
mg 



dx 


(B8) 


Expanding the term in brackets in a power series in {x—Xg) 
jmg, neglecting terms involving the cube and higher powers, 
and noting that 

6 

there is obtained 

(*C fJ-pf 

If equations (6) and (7) are substituted into equation (B9) 
and the indicated integration is performed, ythere is finally 
obtained 


The details of the integration are given in reference 2, page 4. 
The resulting expression is 

^=X An-ry^A„ cos (B5) 


The vertical velocity at a point on the chord of an airfoil 
in cascade is made up of three parts. These are (1) the 
velocitj' v’ induced by its own vortices, (2) the velocity Ai’s 
induced by the vortices of the remaining airfoils of the 
cascade, and (3) the velocity Awj induced by the sources and 
s ink s of the remaining airfoOs. 

By analogy with equation (B5) : 



cos nd 


(B6) 


The vertical velocity induced by the vorticity distribution 
of the other airfoils is given by the expression 


dr' . 

1 r dx Vl-AP ^ 

Vl-APcos= 


(B7) 


where 


^"=-2 Ay)-(2^o'+^4i0 cos ej (BIO) 

(BID 
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It is interesting to note that equation (BlO) can be more 
directly obtained by replacing each of the airfoils of the 
cascade by a single vortex. The important point in such a 
derivation is the proper chordwise location ol this vorte.x at 
the airfoil center of pressure. 

In a manner similar to that in which equation (BlO) was 
obtained it is possible to show that neglecting terms involving 
(x—Xo)(mg to the third and higher powers, 

^=0 (B12) 


In Appendix B of reference 3 the velocity components in- 
duced by a source at the origin are given as follows: The 
velocity component perpendiciflar to the stream Ac, is 


The summation is over the remaining airfoils of the cascade. 
The angle r is measured clockwise between the vertical and 


Q ryi-A/Msin $)~ [ 

2irrl_ 1— A/-sin^<l> J 
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and the induced velocity component parallel to the stream 
Aii, is 

cos $ 


27rr l_Vl— (1— sin^ <t>) 




where Q is the mass flow divided by the free-stream density. 

The vertical velocity induced by the source-sink distribu- 
tions of the remaining airfoils of the. cascade can then be 
expressed as 


Avs=- 


r dQ'Idx r sIl-M- (cos r^) ~| 

7^1 Jo rrn L cos^ Tn J 


and the corresponding horizontal induced velocitj^ is, after 
noting that cos 4>= — sin Tm, 


Au, 


^ m=iJo ^ 


dQ'Idx 


sm T„ 


L 


Vl-^/Tl-M=cos' 


Vm)_ 


dx 


dQ' 


where dx is the source strength over an element of chord 

dx of the cascade airfoil. Substituting for r„ and cos r,„ and 
sin Tm as before, expanding, and neglecting the thhd and 
higher power of (x—xQjmg, there is obtained 


Am, 


A(,1,5=0 

= -5W‘.C *! 


(B13) 


dQ' 

The source-sink distribution is now expressed by the 
following series: 


dQi._ 

dx'' 


--2V{Bd cot 7 T sin nd) (B14) 


On substitution of equation (B14) the integral in equation 
(Bl3) can be evaluated, obtaining, analagously with equa- 
tions (BlO) and (Bl2): 




(Bio) . 


and 


^^==^[{B,’+^B,')-{2B,' + Bd) cose] (B16) 


Since the airfoil must be a closed figure, it follows that 

'^dQ' 


j: 


dx 


dx=Q 


(B17) 


By substitution of equation (Bl4) in equation (Bl7), there is 
obtained 


so that 


2.B/ + B/=0 


^'^■^=~(2Bd+B,') 


(B18) 


V “X®' 


From thin airfoil theory (reference 2) it can be shown that 

+ - (B19) 

where yt is the ordinate of the base jirofile as measured from 
the axis of symmetry in terms of the airfoil chord at the cor- 
responding j: station. 

In the problem of determining the wnill interference in a 
two-dimensional wind tumiel, an evaluation of the influence 
of the image base profiles yields an identical result if the 
limitations of thin airfoil theory are presumed. Lock (dis- 
cussed in reference 3) has evaluated the image base profile 
effects as regards wall interference in the incompressible case 
for profiles not necessarily thin and has found that (in the 
notation of this report) 




;o) 


is a more precise value to replace 2Bd' -[-Bi , where P/, is the 
base-profile pressure coefficient in incompressible flow. 
Values of A computed by equation (B20) for various airfoils 
are given in table I which was taken from reference 3. 
Accordingly, as in reference 3, this value is used in this report 
so that the velocity increment due to the effect of soiu'ccs and 
sinlcs becomes 


Am, Act 

~V~\^ 


(B21) 


It should be noted that properly Lock’s result for incom- 
pressible flow is 

Am, 


IT 


= Act 


(B22) 


However, the expression for T’' is 

V=\\+AU 

so it follows that 


Iq 

V 

Thus for incompressible flow' 
Au 


Am 

F 


iM, / Am 
y=Aa[^l- 


(B23) 


which to the first order in a is still equation (B22) so that 
equation (B21) for compressible flow' follows. From equa- 
tions (Bl2) and (B21) it is clear that the influence at the 
position of one airfoil of the cascade due to all other airfoils 
of the cascade is simply that promoted by the source-sink 
system 

Au_Atr 


-B24) 


Moreover, because Ar, (equation (B13)) is zero, it follow's 
from equation (BlO) 


Ao 

V 


l=-2^^{Aa'+lA/)-{2A,'+Ad) cose] (B25) 



APPENDIX C 


SAMPLE PRESSURE-DISTRIBUTION CALCULATION 

Table II gives all the necessary calculations to arrive at the 
cascade pressure distribution shown in figure 6. It was 
desired to obtain the pressirre distribution for a cascade of 
NACA 4412 airfoils given the values c/=1.0 and ejg=l.QZ 
for comparison with the distribution calculated by the method 
of conformal transformation given in reference 11. 

Since the assumed flow is incompressible, then X, n, and | 
are unity and and ij' are zero. From table I, A is 0.237. 
The parameter <r is 0.220. The hft coefficient for the 
corresponding airfoil in free air is obtained from equation 
(48) as 1.335. By potential theory, for the NACA 4412 


airfoil the free-air angle of attack for Cj= 1.335 is 6.8° and 
is —0.11 so that for the airfoil in cascade, by equation 
(22), ce'=10.4°. 

Values of the pressure coefficient for the isolated airfoil 
were calculated from the free-air values of reference 1 1 using 
the method of reference 5. From these pressures P was 
determined and the P* was obtained from equation (24). 
P* and 1 — P/* (equal to 1 — P/) were then combined by 
equations (43) to give 1 — Pu* and 1 — P^*. The cascade 
pressure distribution was finally obtained using equations 
(44). 


TABLE I.— VALUES OF A FOR VARIOUS BASE PROFILES 
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TABLE II.— CALCULATION OF C.ASCADE PRESSURE DISTRIBUTION OF FIGURE 6 (SEE APPENDIX C) 

[<r=0.220: A =0.2371 
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